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,,^J ' I review three different non-perturbative approaches to the three dimensional 

D , A'^j^-flavor Thirring model: the 1/Afy expansion, Schwinger-Dyson equations, and 

^^ ■ Monte Carlo simulation. Simulation results are presented to support the existence 

of a non-perturbative fixed point at a chiral symmetry breaking phase transition 
for Nf = 2 and 4, but not for Nf = 6. Spectrum calculations for Nj = 2 reveal 
conventional level ordering near the transition. 



1 Introduction 
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CN \ In this talk I will present work done in collaboration with Luigi Del Debbie and 

^^ ■ John Mehegan El. The Thirring model is a field theory of relativistic fermions 

interacting via a contact term between conserved currents. Its Lagrangian is 

written 



C = ^,(^ + m)V'. + l^{i^^l^.^^)\ (1) 

^- 

O . where ipi^ipi are taken to be 4 component spinors with flavor index i running 

from 1 to Nf. Why is this an interesting model to study in three spacetime 
dimensions? Firstly, it appears that the expansion in powers of 1/iV/ is ex- 
actly renormalisable for 2 < d < 4, which in turn suggests an equivalence 
between the strong coupling limit of the Thirring model and the infra-red 

C^ ' behaviour of QEDg. This result is essentially a perturbative one, albeit us- 

ing a different expansion parameter; however non-perturbative studies using 
e.g. the Schwinger-Dyson approach suggest the model may exhibit dynami- 
cal symmetry breaking and have a non-trivial ultra-violet fixed^oint, of the 
type discussed by Miransky in the context of quenched QED4 a, for Nf less 
than some critical Nf,, < 00. At the proposed fixed point there would be no 
small dimensionless parameter in play; this is therefore a natural candidate 
for investigation by Monte Carlo simulation of the lattice-regularised model. 
In the following I will review each of these approaches, before presenting our 
results. Finally, it is worth mentioning that non-trivial fixed-point behaviour 
of QED3 has been proposed to account for non- Fermi liquid behaviour in high 
Tc superconductors S. 

1 



2 Three non-perturbative approaches 

2.1 i/Nf expansion 

The Lagrangian {M may be rewritten using a vector auxiliary field A^: 

£ = ij,{^ + m)^P, + -jl-A^,{l;,-i^i:, + -A^A^. (2) 

In weak coupling perturbation theory, the vector propagator I?^^(fc) is simply 
6^,y\ however to leading order in 1/-/V/ the vacuum polarisation bubble, also 
contributes, with the result that D^,y acquires momentum dependence qEI: 

/ k k \ 1 

Df^Ak) = Um-^ - -^ j i_ii(k2\ + longitudinal, (3) 



with 



^l(m+ -^{k^ ~ 4m2) tan-i (^ 
This form can be examined in two limits. In the UV regime 



The 1/fc asymptotic behaviour has the effect of softening the contact term and 
thus decreasing the superficial degree of divergence associated with higher order 
diagrams containing internal vector lines, such as fermion self-energy, two- 
loop vacuum polarisation, etc. These divergences can be absorbed by a finite 
number of counterterms, essentially due to a novel application of Weinberg's 
theorem, rendering the 1/Nf expansion renormalisableLI. At lowest non-trivial 
order, the corrections coincide with those calculated in the IR limit of QED3LI, 
giving the first hint that the two models may be related. 

In the low-momentum regime, (1^,^) can be used to find a pole condition 
and hence calculate the mass My of the vector bound state. The result depends 
on the dimensionless combination {mg^)0: 

f2TOfl-2expf — ^]] mg2^0: 
Mv={ y^ '^ "^^^ ' (6) 

Thus we identify a smooth crossover from a weakly coupled theory in which 
the vector is a weakly bound // state to a strongly coupled one in which it 



is almost massless. To leading order in l/iV/ the {■ip'yf^ip)'^ interaction is a 
marginal operator. In the strong coupling limit g^ —^ oo a. massless vector 
mediates interaction between conserved currents, suggesting that in this limit 
the theory is identical to QED3. 

2.2 Schwinger-Dyson equations 

The Schwinger-Dyson approach Q is to solve for the full fermion propagator 
Sp^{k) = A{k)i^ + E(fc) self-consistently, by finding integral equations foe 
the functions A and S. The most systematic analysis of the Thirring model H 
exploits a hidden local symmetry of (0) to fix a non-local gauge in which the 
wavefunction renormalisation A{k) = 1. The usual truncations of setting the 
full vector propagator to its leading order 1/Nf form (^) and the full vertex 
rp(p, q) to the bare vertex igjfj,/ y^Nf (i.e. the ladder approximation) are also 
employed. The resulting SD equation for I](fc) can be solved analytically in 
the limit g^ —>■ 00. 

In the chiral limit m — > 0, a non-trivial solution for S(A;) is found for 
small values of Nf, signalled by a non- vanishing order parameter, the chiral 
condensate (ipip)'- 

{U) ex A^ exp(-37r/^^ - 1), (7) 

with A a UV cutoff. This solution exists for Nf < Nfc = US/Stt"^ ~ 4.3, 
implying that in the strong coupling limit, a separation of IR and UV scales 
allowing a continuum limit can only be made as TV/ ^p- Nfc- The essentially 
singular behaviour is reminiscent of quenched QED4I3. For finite g^ there is 
no analytic solution, but it has been argued that the critical curve Nfdg) is 
smoothB, implying that non-trivial UV fixed points as g^ — > g1{Nf) may exist 
for integer Nf < Nfc- The Thirring interaction has thus become relevant. It 
should be noted, however, that a different truncation procedure yields Nfc — 
00 13. It is therefore desirable to try an alternative non-perturbative approach 
to see if fixed points can be identified, and to get an independent estimate of 
Nfc 

2.3 Lattice formulation 

The lattice action we have used is 



S 



^ X! ^' (^^'^'^ (^) [^* {x + i^) - Xt{x - fi)] + 



2 

„2 

m^Xt{x)Xi{x) + ^ X! Mx)Xi{x + fi)x3{x + fi)Xj{x) (8) 

xi xfiij 
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Figure 1: 



Figure 2: 



where i runs over N flavors of one-component staggered ferniion fields XiXj 
and r]fj,{x) are the space-dependent Kawamoto-Smit phases which are lattice 
analogues of the 7-niatrices. We have simulated the models with N = 1,2 and 
3, proceeding by introducing a vector auxiliary on lattice links by analogy with 
(0). Alternative formulations are possibleliJ; we have chosen (g) so that there 
is no gauge symmetry, and the highest interaction is four-point. 

There are two aspects of the formulation (0) which require further discus- 
sion. Firstly, and slightly unexpectedly, the lattice interaction current is not 
exactly conserved, because it couples to iAf^ and not exp(i^^). This means 
that the lattice vacuum polarisation tensor is not transverse, which at leading 
order ia 1/-^/ results in an additive renormalisation of the inverse coupling 



f: 



9l = 



1- 



h^' 



(9) 



9l 



We should thus expect to find the strong coupling limit at 

1/g'^ -^ 2/3 (for 1/g^ < 2/3 the model is not unitary). It is clear that for small 

Nf, identifying the limit g^ —> 00 numerically will be difficult in practice. 

Secondly, when the lattice action (H) is recast in terms jXif physical fields 
q, q which have explicit spin and flavor indices, then we find El 
(i) Each flavor of x corresponds to 2 physical flavors, via fermion doubling, i.e. 



Nf = 2N. 

(ii) In addition to the expected interaction term {qi{"ff_i (X 
second 2x2 matrix in the product acts on flavor indices, 



(10) 

l)(7i)^, where the 
) contains other 



Table 1: Fits to the Equation of State 



_^ TV/ = 2 Nf ^4 

'^Jgl 1.92(2) 0.66(1) 

5 2.75(9) 3.43(9) 



interactions such as {(jii'jzj^ <E) l)?;)^ ^^id even non-covariant ones such as 
(9i(747j/7Ai'8)T*)gi)^. These undesired terms are not suppressed in the hmit lat- 
tice spacing a ^ 0. Since a strongly-coupled theory requires a non-perturbative 
regularisation, we must take a robust attitude to this problem, and determine 
the fixed point behaviour, if any, by numerical measurement of critical scaling 
and spectral quantities. 

3 Numerical results 

3.1 Critical scaling 

We have Monte Carlo data fronp. lattice sizes 8^, 12^^ and 16^, and a range of 
bare masses ma — 0.05, . . . , O.OlEl. To determine the critical properties requires 
an extrapolation to the chiral limit m ^ 0; this is best done globally, so in the 
neighbourhood of the chiral phase transition we assume an equation of state 
of the form 

m = A(^-\]{i,^)^-^/'' + B{i,^/jy. (11) 

Here S and /3 are standard critical exponents, which in a mean field treatment 
would assume the values 6 = 3, (3 = ^. It follows that in the mean- field case a 
plot of (ipip)'^ vs. m/ {ijjip), known as a Fisher plot, would yield straight lines 
for constant l/g^ which would pass through the origin at the critical coupling, 
and above or below the origin in the super- or sub-critical region respectively. 
In Figs. 1 and 2 we show Fisher plots for Nf = 2 in the neighbourhood of 
l/gl ~ 1.9 and Nj = A near 1/g^ ~ 0.6. The hues denote our best fit|±o ([u]) 
assuming S — 1/(3 = 1, which follows from the ladder approximation 113. The 
curvature indicates a slight departure from mean-field scaling. Our fits, which 
include a finite volume scaling analysis for Nf — 2, are presented in Tab. || 

The remaining exponents can be derived using the constraint between S 
and /3, and hyperscaling. The most significant result is not so much that the 
value of S departs from the mean-field value in each case, but rather that the 
values from each model are distinct. This is consistent with the Nf = 2,4 
theories belonging to different universality classes. 
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A Fisher plot for Nf — 6, shown in Fig. 3, shows no evidence for a chirally 
broken phase, and no fit of the form ( pj] ) is possible. Thus the data, from 
admittedly small lattices to date, strongly suggest 



Nf, < 6. 



(12) 



3.2 Spectroscopy 

We have measured the physical fermion and bound state masses for Nf = 2 on 
a 16'^ lattice with ma = 0.01. Fig. 4 shows fermion, pion and scalar masses vs. 
1/5^, showing a sharp crossover near 1/g^ ~ 1.9. In the symmetric phase both 
pion and scalar look like weakly bound // states, while in the broken phase 
the pion is the lightest particle, consistent with Goldstone's theorem. There 
is no evidence for the absence of light bound states in the symrrietric phase, 
predicted if the fixed point is described by an essential singularity t3. We have 
also examined the vector channel, where the simplest interpolating operator 
projects onto states with both vector {<pi^q) and axial vector {q^^^^q) quantum 
numbers. Surprisingly, we found both states to be light in the symmetric phase, 
even though to leading order in 1/iV/ there is no bound state in the axial vector 
channel. This is a hint that the l/A^/ expansion may be failing to describe the 
physics even of the symmetric phase near the fixed point. Measurements in the 
broken phase are much noisier, but there is tentative evidence that the vector 
mass increases sharply. 



4 Summary 

Our main results are that for Nf — 2,4 there is clear evidence for a continuous 
chiral phase transition, with scaling described by distinct critical exponents 
in each case. No such transition is seen for Nf = 6, implying Nfc < 6. For 
Nf = 2, we see conventional level ordering of the spectrum in fermion, pion 
and scalar channels across the transition, and light states in the symmetric 
phase in both vector and axial vector channels. 
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